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Abstract. 

During the past fifteen years, inhomogeneous cosmological models have been put 
forward to explain the observed dimming of the SNIa luminosity without resorting 
to dark energy. The simplest models are the spherically symmetric Lemaitre-Tolman 
(LT) solutions with a central observer. Their use must be considered as a mere first 
step towards more sophisticated models. Spherical symmetry is but a mathematical 
simplification and one must consider spherical symmetric models as exhibiting an 
energy density smoothed out over angles around us. However, they have been taken 
at face value by some authors who tried to use them for either irrelevant purposes or 
to put them to the test as if they were robust models of our Universe. We wish to 
clarify how these models must be used in cosmology. We first use the results obtained 
by Iguchi and collaborators to derive the density profiles of the pure growing and 
decaying mode LT models. We then discuss the relevance of the different test proposals 
in the light of the interpretation given above. We show that decaying-mode (parabolic) 
LT models always exhibit an overdensity near their centre and growing-mode (elliptic 
or hyperbolic) LT models, a void. This is at variance with some statements in the 
literature. We dismiss all previous proposals merely designed to test the spherical 
symmetry of the LT models, and we agree that the value of Hq and the measurement 
of the redshift drift are valid tests of the models. However, we suspect that this last 
test, which is the best in principle, will be more complicated to implement than usually 
claimed. 
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1. Introduction 

The dimming of the SN la luminosity discovered by Riess et al. (1998) and Perlmutter 
et al. (1999) was first interpreted in a standard FLRW cosmology as due to the effect of 
a cosmological constant or a dark energy component. Since then, a number of alternate 
solutions have been proposed. Among them, exact inhomogeneous models with no dark 
energy component were put forward shortly after the release of the first supernova results 
(Dabrowski and Hendry 1998; Pascual-Sanchez 1999; Celerier 2000; Tomita 2000, 2001; 
Iguchi et al. 2002). Then, after some neglect, they have started to be reconsidered 
(Rasanen 2004; Celerier 2007). 

For simplicity sake, most of the first models exhibited spherical symmetry around 
the observer. This was actually the property used to give the first general analytical 
derivation of the way dark energy could in principle be explained away by local 
inhomogeneities (see Section 2 of Celerier 2000). However, we wish to dismiss once more 
a biased interpretation that has been too widely spread since then in the cosmological 
community. In contrast to what has been advocated by the author of the present paper 
in the above cited article and by others afterwards, the use of this kind of models should 
not imply that the observer is physically located at the centre of a spherically symmetric 
universe. At the time these models were first proposed, the available cosmological data 
seemed of course to point to an isotropic universe around us on all cosmological scales. 
However, since then, more and more precise observations have shown structures on larger 
and larger scales and we can now conclude that the only nearly isotropic feature in the 
Universe is the CMB. Close to the last scattering surface, radiative pressure cannot be 
neglected and dust models such as LT models are no longer relevant. 

We therefore propose that one should consider these inhomogeneous spherically 
symmetric models as exhibiting an energy density smoothed out over angles around 
us (Celerier et al. 2010; Bolejko and Sussman 2011; Bolejko et al. 2011), i.e. only 
the radial inhomogeneities are taken into account. This was exemplified by Bolejko 
and Sussman (2011), who showed, with a realistic Szekeres model, that spherical void 
structures approximate configurations that can emerge by coarse-graining and averaging 
a large-scale region of a lumpy universe in which the density distribution is far from 
spherical. 

One should consider that going from the everywhere homogeneous FLRW model to 
the spherically symmetric LT model merely amounts to switch the use of a solution with 
three Killing vectors (FLRW) for one with only two Killing vectors (LT). One symmetry 
(the radial one) is deleted. This was a mere first step to dealing with the cosmological 
constant problem. 

Spherical symmetry around the observer is, of course, a mathematical simplification 
that is disregarded by more complex models such as Swiss-cheese models with spherical 
holes (Brouzakis et al. 2007; Marra et al. 2007; Biswas and Notari 2008; Brouzakis et 
al. 2008; Marra et al. 2008), meatball models (Kainulainen and Marra 2009), Szekeres 
Swiss-cheese models (Bolejko and Celerier 2010) etc. Thus, these models must not be 
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taken at face value and some of the proposals designed to test them as if they were exact 
representations of our Universe, are irrelevant. 

Another strange argument against these kinds of models is that they are accused of 
being non-Copernican, while the standard spatially homogeneous models are credited 
to obey the so-called Copernican principle. However, if one considers an FLRW model 
from a fully relativistic point of view, i.e. if the location of the observer is considered not 
only in space but rather in the four- dimensional spacetime, the model can no more be 
called Copernican. In reality, the location of the observer in space is not special, but his 
temporal location in the expansion history of the Universe is. In the framework of the 
dark energy paradigm, this issue is known as the coincidence problem. Therefore, when 
such an argument is put forward to suggest that FLRW models should be preferred to 
spherically symmetric ones, this seems rather inconsistent, even from a philosophical 
point of view. 

The spherically symmetric models most commonly used to attempt to explain away 
dark energy are those of the Lemaitre (1933) - Tolman (1934) (LT) class. These are the 
ones we discuss in this article. 

In Sec. El we recall the equations defining the LT solutions, mainly to set the 
notations that are used here. In Sec. [3l the shape of the mass density profiles of 
different classes of LT models are studied and some new clarifications are given. Section 
mis devoted to the discussion of different tests of these models proposed in the literature 
and, in Sec. [5l we present our conclusions. 



2. The Lemaitre — Tolman model 



We recall the equations needed in the following to study the main features of the LT 
models. We adopt widely used notation defined in, e.g., Krasinski (1997), Plebanski 
and Krasinski (2006) or Bolejko et al. (2010). 

The LT model is a spherically symmetric non-static solution of Einstein's equations 
with a dust source. Its metric, in comoving coordinates, synchronous time gauge and in 
units in which c = 1, is 

ds^ = dt" - — ^^^dr^ - R^{t,r){M^ + sin^ ^9 df/^^)^ (1) 
1 + 2hi(r) 

where E{r) is an arbitrary function of integration determining the curvature of space at 
each r value, a prime denotes a derivative with respect to r and R{t, r) obeys, in units 
in which G = 1, 

•9 2M A , 

R^ = 2E + — + -R\ (2) 

where a dot denotes a derivative with respect to t and A is the cosmological constant. 
Equation (j2]) is a first integral of one of Einstein's equations, and M = M{r) is another 
arbitrary function of integration representing the gravitational mass contained within 
the comoving spherical shell at a given r. The mass density is 
2M' SnG 
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In the following, we assume A = 0. In this case, Eq. 
and the solutions are: 



can be solved explicitly. 



when E <{) (elliptic evolution) 

M 



R{t,r) 



rj — sm rj 



:-2E) 



M 



1 — cos rj) , 

(i-^B(r)), 



(4) 
(5) 



where ri{t,r) is a parameter; 
when E = (parabolic evolution). 



R{t,r) 



1/3 



when E > (hyperbolic evolution) 

M 



R{t,r) 



2E 



(cosh 7] — 1) 



sinh f] — rj 



{2E) 



3/2 



M 



-{t-tB{r)). 



(6) 



(7) 



(8) 



All the formulae presented so far are covariant under arbitrary coordinate 
transformations f = g{r), such that r can be chosen at will. This means that one 
of the three functions -E'(r), M(r), and ^^(r) can be fixed at one's convenience by an 
appropriate choice of g. 



3. Shape of the mass density profiles of LT models reproducing the ACDM 
luminosity distance-redshift relation 

Most of the density profiles adopted in the literature to construct LT models of a universe 
with no dark energy that are able to reproduce a significant part of the cosmological 
data are of the void type (see, e.g., for reviews, Rasanen 2006, Celerier 2007, Bolejko et 
al. 2011, Marra and Notari 2011). It has been proven, however, that a void model is 
not mandatory to reproduce the main observables of the ACDM model (Celerier et al. 
2010). 

We suspect that the choice of a void is dictated by the reproduction of an accelerated 
expansion being easier to grasp within such a model. To provide indeed an intuitive 
understanding of how inhomogeneities can mimic an accelerated expansion, many 
authors have cited the very pedagogical model proposed by Tomita (2000, 2001). This 
model is composed of a low-density inner homogeneous region connected at some redshift 
to an outer homogeneous region of higher density. Both regions decelerate, but, since 
the void expands faster than the outer region, an observer located inside this void 
experiences an apparent acceleration when she looks at a source located in the outer 
region. 
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There is some confusion in the hterature about which models are voids. For 
instance, the LT example of Celerier (2000) and the E = {] LT case of Iguchi et al. 
(2002) have often been presented as void models, which they are not as we show below. 

We have seen above that the r coordinate can be chosen at will by fixing the form 
of one of the three first integrals of Einstein's equations, which are functions of r. Here, 
we make the rather common choice of M(r) = M^r^ . This choice also has the advantage 
that it is the one used in the above cited articles, namely Celerier (2000) and Iguchi et 
al. (2002), whose interpretations we intend to specify correctly here. 

3.1. The parabolic (flat) E = case 

We first consider the parabolic (flat) case with E = 0. With the gauge choice described 
above, Eq. ([6]) becomes 

Inserting into Eq. [3]) this expression, its derivative, and the derivative of M(r) with 
respect to r, we obtain 

_ 4 

~ {t-tB){3t-3tB-2rt'sy 
We then calculate the derivative of p with respect to r 



np' (t - ts) (St's + 2rt%) - 2rt 



4 {t-tBy{3t-3tB -2rt'^ 



B) 



\2 ■ 



To see whether the observer, with time coordinate t^, is located on a density hill or 
in a void, we take the limit of the above expression when r — )■ 0. With our gauge 
choice, we can write tB{r) = TiT + r2r^ + 0{r^), where we have renormalized by 
translation to be tsir = 0) = 0. Therefore, near the centre, we have {to — tB) — to, 
(3to -3tB- 2rt'^) -> 3to and 2rtot'^ - 2rt% < Sto^B StoTi. Thus, we obtain 

32t 

/€p'(t = to,r->0)^— 3^. (12) 
The same limits applied to Eq. ffTOl) give 

Kp{t = t^,r ^Q)^ ^. (13) 

This implies that, at the observer, the density profile is finite and exhibits a cusp. 

Since we can, with no loss of generality, choose to > 0, the direction of its slope at 
the centre, which determines whether the observer is located in a void or on a density 
hill, is given by the sign of Ti. 

We can first note that, to avoid shell-crossing, ri < (Bolejko et al. 2010). However, 
another check can be completed. 

The sign of ri is the same as that of t'^ when r — )■ 0. Iguchi et al. (2002) studied the 
two pure classes of A = LT models which reproduce the luminosity distance-redshift 
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relation of the ACDM model, one with = and the other with ts = 0. One can see 
in their Fig.2 that, for E = 0, <\tB/(\z{z = 0) < 0. 
We now use 

t' = '^ = (14) 
^ dr dz dr 

For a radial incoming ray (Bondi 1947), 

dz _ (1 + z)R' 
dr ~ VI + 2^' 
In the case = 0, 
dz 
dr 

with 



(15) 



l + z)E!, (16) 



^, '4M,y/=3«-3<^ (17) 



At the observer {t = t^.r ^ Q), 



K,^['-^T>^- (18) 



3to , 

We thus see that the sign of t'^ near the observer, which is that of Ti, is the same as 
that of (itB/(\.z{z = 0), i.e., a minus sign. 

This implies, therefore, that the observer is located on a density hill and not in a 
void. As previously noted by Celerier et al. (2010), this feature of the density profile 
only exists on the t = to spacelike hypersurface and is thus not directly observable (see 
however the discussion about the redshift drift in Sec. Hj). 

We therefore wish to stress that the example studied by Celerier (2000) and the 
first case considered by Iguchi et al. (2002) are not void models but overdense ones, in 
contrast to some incorrect statements in the literature. 



3.2. The simultaneous bang ts = case 

We now consider the simultaneous bang case with ts = 0. It is well-known that, with 
only one degree of freedom left when choosing two of the arbitrary functions of r defining 
the model, M and ts, the supernova data completely fix the features of the last function 
E{r). However, we subsequently see that the sign of the derivative of the mass density 
near the observer is the same in both cases, E < and E > 0. We therefore study 
separately the two cases below, without bothering about whether we know the sign of 
E{r) at the observer, which is not our purpose here. 
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3.2.1. The elliptic E < case We first calculate the energy density p{r) at the observer. 
For this, we need the expression of R'. Differentiating Eqs. (jlj) and ([5]) with respect to 
r and setting = and M = Mqt^, we obtain 



2 r 



R' 



{=2E) 
3t 

~ (-2^)1/2 

Inserting this expression and 

6 _ Mor^ 
up 



3 + 



2rE' 



E' + 



-2E) 
-2E) 



1 — cos rf) 



1 — cos 7] 



(19) 



in ([3]), the energy density is given by 



{-2Ef 
X (1 — cos rj) — 



(1 — cos 77)^ 
3t 



Mor^ 



3+ 



;-2E)V2 



E' + 



-2E) 



2rE' 

h2E) 
sin?7 



1 



cos rj 



(20) 



We can always, without lack of generality, since r/ is a function of t and r, renormalize 
to at the observer such that, at r = 0, 1 — cos rj and sin 77 should be nonzero. In this 
case, since E = Eqt^ + 0{r^) and thus E' = 2Eor + 0{r^), when r — )■ 0, p exhibits a 
finite value. 

We now wish to study the shape of the density function near the observer. To 
complete the task, we cannot use the same method as in the E = case, since this 
would not allow us to draw any definite conclusions. Therefore, we adopt the following 
reasoning. 

From Fig.5 of Iguchi et al. (2002), we see that dp/dz > at the observer. Now, 
, dp dz 
^ dz dr ' 
and for a radial incoming ray, 

dz + z)R' 
dr 



(21) 



(22) 



^/l + 2E' 

For r — )■ 0, we have z — )■ 0, 1 + 2E — )■ 1, and therefore dz/dr — )■ R'q. To calculate i?', we 
differentiate the expression of R' given by Eq. (fT9l) with respect to t. This gives 



R' 



X 



-2E) 



1/2 



E' + 



r 

:-2E) 



3 + 



2rE' 

sin rj 
1 — cos rj 



+ 



(-2^)1/2 
3t{-2E) 



Mor3 



E' + 



i-2E) 



'1 — COST]) 



2 ■ 



(23) 



To compute R'q at the observer, it is sufficient to keep the first term in the expansion of 
E{r) in powers of r. We thus obtain 



-2E,Y/'- 



sm rjQ 



1 — cos rjo 

Moreover, from Eq. ([5]), we can write near the centre 

-2Eof/\ 



r]o - sm r]o 



Mn 



-tn. 



(24) 



(25) 
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Since (— 2£'o), to ctud Mq are positive, the above equation implies that rj^ > sinT^o, and 
therefore < r^o < which gives sinr/o > 0. Since we have renormahzed to so that 
1 — COSTCO 7^ 0, we have 1 — cos?7o > 0. From Eq. ( l24l) . we thus obtain i?Q > and 
therefore, (d2;/dr)o > 0. Hence, Eq. (12T!) imphes that Pq has the same sign as (dp/d^)o, 
which is strictly positive as Rq. In this case, the observer is thus located in a local void. 

For an elliptic evolution, we can therefore conclude that the requirement that the 
bang function should be a constant (always renormalizable to zero), which has been 
widely assumed, as in, e.g., Yoo et al. (2008), Zibin (2008) and Zibin et al. (2008), 
implies necessarily a void. 

This also demonstrates that pg is non-zero (p is non-differentiable at the centre). 
This was previously noted by Clifton et al. (2008). In all cases, we claim that this 
property cannot be used as an argument to dismiss the model, since other features of 
this kind can be found in nature, e.g., on the surface of the Earth (Bolejko et al. 2011). 

3.2.2. The hyperbolic E > case We again first calculate the energy density p(r) at 
the observer. Differentiating with respect to r Eqs. ([7]) and ([8]) where we have set = 
and M = Mqt^, we obtain the expression of R' 



, Mor^ ( rE'\ , , 



3t / 2E\ sinhr] 

+ (2E)i/2 V V; coshr^- r ^ ^ 



Inserting this expression and Eq. ([7]) into Eq. ([3]), the energy density is given by 

(coshr] — 1) 



^ — ° :(cosh r] — 



2 



Kp {2E) 

^ 3t 2E\ sinhi] 







. 2^ 





(27) 

(2E)V2 \^ r J coshr] - 1_ ' ^ ^ 

As in the previous subsection, we can renormalize to at the observer such that, at 
{t = tQ,r = 0}, cosh?7o — 1 and sinh?7o should be nonzero. Therefore, p exhibits a finite 
value there. 

We now study the sign of p'q = p(t = to,r — )■ 0). We have already noted that, 
according to Fig.5 of Iguchi et al. (2002), dp/dz > near the observer. Using Eq. (122|) . 
we have shown above that for {t = to,r ^ 0}, dz/dr — )■ i?Q. The calculation of R' here 
gives 



R' 



{2Efl^ f^_r_E;\^ 3 (j^, _ 2E\ 



r V E J (2E)V2 \^ r J 

=i„h„ BE 1 (28) 



coshi] — 1 Mor^ \ r J (coshr/ — 1)^ 

To compute i?Q, we only keep the first term in the expansion of E{r) in powers of r. We 
thus obtain 

p/ _ /op u/2 sinhr/o 

coshr^o-1- ^ ^ 
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From Eq. ([S]), we can write 

smh r]o-r]o = — — to- (30) 

Mq 

As in the elliptic case, the right-hand side of this equation is positive, which implies that 
sinhr^o > ''70; hence r^o > and sinhr^o > 0. Since coshr^o — 1 > 0, we obtain Rq > 0. 
From (12T|) . we see that pg = (dp/dz)o-Ro > 0- ^^^^ case, the observer is thus also 
located in a local void, which exhibits a cusp at the origin. 

As stressed by Krasihski et al. (2010), since an L-T model is fully specified by two 
physical functions, a random combination such as£' = 0ort5 = and an arbitrarily 
chosen Dl{z), such as the one prevailing in the ACDM picture, may well produce an 
LT model with unrealistic features. 

Since the argument that decaying modes must be avoided to construct a model 
compatible with the CMB requirements remains disputed (Yoo et al. 2008; Zibin 
2008; Zibin et al. 2008; Celerier et al. 2010), one should consider that intermediate 
models with both growing and decaying modes can also be used to solve the dark 
energy problem. This was the result obtained by Celerier et al. (2010), where the 
arbitrary functions determining the LT models capable of reproducing a few features of 
the ACDM model were computed without any prior assumptions about their form. In 
these LT models, the shape of the current density profile is not a void, but a hump. 

We also recall here that the function p(r) considered above is not directly observable. 
It exists in the space t = t^ oi events simultaneous with our current instant in the 
cosmological synchronization, i.e. it is in a space-like relation to us. This is also the 
case for the voids studied by Alnes et al. (2006), Garcia-Bellido and Ilaugb0lle (2008a), 
and Yoo et al. (2008). 



4. Test proposals and ruling out claims 

We discuss the main tests and invalidate some claims that have been put forward in 
the literature. We do not intend to give a review of the subject (this can be found, 
e.g., in a recent special focus issue of Classical and Quantum Gravity on inhomogeneous 
cosmological models and backreaction and averaging in cosmology (Anderson and Coley 
2011), but only to set the record right about the newest proposals and claims which we 
find relevant to illustrate our purpose. 

To begin with, we wish to stress that all the published tests have been applied to 
very few particular models of the infinitely larger class of inhomogeneous spherically 
symmetric models and that, therefore, claims that some test or other rules out the 
most favoured void model is premature. Moreover, non-void models have seldom been 
considered but seem to pass some of the tests more successfully. 

However, the main drawback of most of these proposals is that they are designed 
to test the spherical symmetry of the models as if it were an actual physical property of 
the Universe and not a mere mathematical simplification as we have explained in Sec. 

m 
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4.1. The kinematic Sunyaev-Zel'dovich effect 

The use of the kinematic Sunyaev-Zel'dovich (kSZ) effect (Sunyaev and Zel'dovich 1972) 
to test the spherical symmetry of the Universe around the observer was first proposed 
by Garci'a-BeUido and Haugb0lle (2008b). This proposal involves trying to go beyond 
the observer's light cone to place constraints on the models. The idea is that, if we 
are located very near to the centre of a spherical void, we can observe distant sources 
that are off-centre, which therefore corresponds to observing a large dipole in the CMB 
spectrum. Such a dipole would manifest itself to us through the kSZ effect. By using 
available measurements for nine galaxy clusters, the authors claimed that they were able 
to put constraints on a particular void model. 

In recent articles (Zhang and Stebbins 2011, Zibin and Moss 2011, Bull et al. 2012), 
the kSZ effect has been once more applied to test different particular void models. 

Besides the many loopholes described by Zibin and Moss (2011), the main drawback 
of the use of the kSZ effect is that, since we are not at the centre of a spherically 
symmetric universe, the observer in the distant source is not off- centre. Therefore, 
she is not supposed to observe a large dipole in the CMB and no kSZ effect is therefore 
observed, in agreement with the measurements. This point was stressed in the discussion 
of Bull et al. (2012). 

Conversely, if an observer in a distant cluster or free electron cloud were also to 
use an LT model to analyse her cosmological data, she would find herself apparently 
"at the centre" of a void or a hump. If she were to apply the same reasoning as that 
given above for the use of the kSZ effect, she would assume that we are far off-centre in 
her spherically symmetric Universe and that we should see a large dipole in our CMB 
spectrum. However, this is not the case and the value of the dipole that we measure has 
allowed us to place stringent bounds on the departures of our location from the centre 
of a LT model (some tens of Mpc), which are much smaller than those of the remote 
clusters or free electrons used by Garcia-Bellido and Haugb0lle (2008b) and Zhang and 
Stebbins (2011). We therefore see that the above reasoning, which cannot be applied 
both ways, is unsound. 

4-2. Cosmic parallax 

The cosmic parallax effect was proposed by Quercellini et al. (2009) and Quartin and 
Amendola (2010) to put bounds on the departure of a possible off-centre observer from 
the symmetry centre of LT models. This effect is based on, in these models, the off- 
centre observers seeing an anisotropic space. If the expansion is anisotropic, the angular 
separation between two distant sources varies in time, thereby inducing a cosmic parallax 
effect, analogous to the well-known stellar parallax, save that this cosmic parallax is 
caused by a differential cosmic expansion rather than an observer's movement. This 
effect should be measurable by future space missions such as GAIA or SIM. 

However, if this effect can provide an interesting measurement of the general 
anisotropy of our Universe, it is not adapted to say anything about a universe 
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approximated by LT models. As explained in Sec. [H these models are only applicable 
to a central observer, since they are isotropic hy construction. 

4-3. Spectral distortion of the CMB 

Constraints on void models using the spectral distortion of the CMB have been proposed 
in the literature (Caldwell and Stebbins 2008; Moss et al. 2011; Zhang and Stebbins 
2011; Zibin and Moss 2011). The philosophy is the same as that underlying the proposal 
to test LT models with the kSZ effect. A Compton scatterer at a given redshift on the 
observer past light cone observes an anisotropic CMB that is reflected back to us in 
the form of spectral distortions, i.e. deviations from a black body spectrum. This 
y distortion can be compared to the measurements by FIRAS on board the COBE 
satellite (Fixsen et al. 1996), which implies that y < 15 x 10^^. Tight constraints were 
found by Caldwell and Stebbins (2008), who however did not use an LT model but a 
Hubble bubble one. Other studies completed with pure growing mode LT models found 
much weaker constraints (Moss et al. 2011, Zibin and Moss 2011). 

Zibin (2011) analysed a specific LT model with pure decaying mode, i.e. a local 
hump as shown in Sec. 13.11 A generic feature of LT decaying modes is the presence of 
blueshifts. This means that, in an exact LT universe, the scatterer whose last scattering 
surface intersects a decaying mode would observe the CMB temperature in the direction 
of the decaying mode to be greater than the actual temperature at recombination. 
Therefore, it would see a strongly anisotropic CMB and produce significant Compton 
spectral distortions. A conclusion of Zibin (2011) is that pure decaying modes that are 
significant today and wider than 0.02 in redshift are very likely ruled out by this test. 

However, this test exhibits the same loophole as those we previously examined, i.e. 
it is based on the assumption that the Universe is exactly LT. This hypothesis is not 
needed to validate the use of such models as mere approximations, as we wish to do 
convince here the reader. Therefore, these three last tests only tell us what we already 
know from observations, that the Universe is not spatially spherically symmetric around 
us. 

4.4- The value of Hq 

In the context of LT models, the most stringent constraint on the present value of 
the Hubble parameter, Hq, comes from the fitting of the CMB power spectrum. This 
spectrum is the result of, first, the imprint of the primordial perturbations onto the last 
scattering surface, and second, the geometry of the Universe between this surface and 
ourselves, which bends the photon trajectory and relates length scales at last scattering 
to angles measured on our sky. Since perturbation theory applied to an homogeneous 
FLRW solution is a good approximation for the Universe at last scattering, the simplest 
solution, usually found in the literature, is an inhomogeneous model where an inner 
LT region is matched at some redshift to an Einstein-de Sitter (EdS) background. 
Here, the features of the CMB spectrum are determined, first, by the EdS cosmological 
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parameters, the expansion rate up to last scattering and the relative components of 
the Universe, second, by the geometry of the model, which differs locally from the 
standard one owing to the properties of the LT model used. Since we do not have a 
well-developed perturbation scheme for LT models, the simplified method used at the 
dawn of inhomogeneous cosmology using these models was as follows: first, calibrate 
the model by fitting the first peak, then, compute the rest of the spectrum as usual with 
one's favorite CMB code applied to the EdS model. That it is widely beheved that a 
LT model reproducing the CMB spectrum exhibits a value of Hq too low with respect 
to the last measured one (Riess et al. 2011) comes therefore from the following wrong 
reasoning, based on this simplified scheme. 

In the standard ACDM model, if the baryon fraction flh is fixed, e.g., by 
nucleosynthesis, there is a strong degeneracy between the fraction of total matter 
(baryon + dark matter) Qm and the Hubble parameter at the observer Hq for a given size 
of the angle 6 subtended today on the CMB sky by the sound horizon at recombination. 
If an LT model is matched to an EdS model at some redshift much lower than that of the 
last scattering surface, we find that — 1 for the CMB. Thus, the size of 9 depends 
essentially on Hq. Now, the ratio of the comoving distance to the recombination dc 
to 9 is linked to the scale of the CMB first peak, which is precisely known from the 
measurements derived from WMAP data. In a given LT model, dc and 9 are easy to 
compute. Therefore, the value of Hq depends on dc{z) = (1 + z)R[t{z),r{z)], i.e., on 
the two arbitrary functions of the model defining R{t,r). If the background is assumed 
to be EdS and the bang time function of the LT model is assumed to be constant, then 
in order to fit a good shape to the CMB power spectrum a low value of the expansion 
rate is required. This is because the proper shape of the power spectrum demands that 
Qj^h'^ Ki 0.13. Therefore, if one assumes that Qm — 1, then one gets h fa 0.4. To fit the 
supernova data, one needs a fluctuation of the expansion rate of amplitude 5h ~ 0.1—0.2 
(Enqvist and Mattsson 2007, Bolejko and Wyithc 2009), so this implies that the local 
expansion rate is low, e.g. Hq ^ 45 km s~^ Mpc^^ (Moss et al. 2011). 

However, it has been shown by Clifton et al. (2009) and Bull et al. (2012) that the 
local Hubble rate is sensitive to the local bang time and that, by making our region of 
the Universe younger we can increase Hq up to the measured value. 

Another way of increasing the value of Hq was proposed by Biswas et al. (2010), 
which is the inclusion of a nonzero overall curvature and a variation in the density profile 
of the void. Here, the best-fit displayed profiles are compatible with ifo = 62.3 ±6.3 km 
Mpc~^, which is the HST value of Sandage et al. (2006). However, the higher value 
of Hq measured at the time, equal to 74.2 ±3.6 km s""*" Mpc~^ (Riess et al. 2009), could 
only be obtained at the price of reproducing less accurately the other observations. 

Some authors have tried more recently to go deeper into the CMB issue in void 
models. Clarkson and Regis (2011) demonstrated that void LT models fitting the 
supernova data and exhibiting a high Hq (= 70 ± 10 km s~^ Mpc~^) can reproduce the 
detailed TT and EE angular power spectra for a flat ACDM model up to 1=2500. This is 
obtained while allowing for the dynamical effect of radiation but ignoring the Integrated 
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Sachs- Wolfe effect. Nadathur and Sarkar (2011) considered a primordial power spectrum 
that was not assumed to be scale-invariant. They were thus able to simultaneously fit 
the SNIa data as well as the full CMB power spectrum, while satisfying constraints from 
local Hubble measurements, primordial nucleosynthesis, and BAO. We note, however, 
that the value for Hq used here is that obtained by the Hubble Key project i/o = 72 ± 8 
km s~^ Mpc~^ which is a little smaller than the latest published value of Hq reproduced 
below. 

However, Riess et al. (2011) claimed that this latest measurement of Hq — 73.8±2.4 
km s~^ Mpc~^ that they made with the Hubble Space Telescope (HST) and the Wide 
Field Camera 3 (WFC3) was likely to provide a strong rebuff to giant void models with 
a central observer. As we have seen above, this is only true for some LT models, but 
not quite all of them, e.g., those of Clarkson and Regis (2011) and Nadathur and Sarkar 
(2011). 

Finally, Romano (2011), using a local redshift expansion for the luminosity distance 
and a constraint on the age of the Universe, showed that the parameters defining a 
general LT model give them enough freedom to enable them to agree with any value 
of Hq. Even if their analysis seems to lack sufficient data to completely constraint the 
models, and in particular constraints from the CMB, this agrees with the results found 
by Celerier et al. (2010), where LT models reproducing a few features of the ACDM 
model were computed with no a priori assumption about the form of the density profile. 
Two separate sets of data were used - the angular diameter distance together with 
the rcdshift-space mass density and the angular diameter distance together with the 
expansion rate - both defined on the past null cone as functions of the redshift. The 
result is not a void but a hump profile and the measured value of Hq is obtained in both 
cases. 

We therefore claim that the value of Hq measured by Riess et al. (2011) only rules 
out some LT universe models but not all of them, be they a void or an overdensity. 

4.5. Redshift drift 

In our opinion, this is the most robust method for discriminating between different 
models reproducing the supernova data, as shown for spherically symmetric and 
homogeneous models by Uzan et al. (2008). 

The redshift drift is the temporal variation in the redshift of distant sources when 
the observation of the same sources is done at different observer's proper times in an 
expanding universe. This was first considered by Sandage (1962). In FLRW models, 
when the expansion of the universe decelerates, all redshifts decrease with time. In 
models where the expansion is recently accelerating, as in the ACDM model, sources 
with redshifts < 2 exhibit a positive redshift drift. Uzan et al. (2008) proposed to use 
this effect to test the so-called Copernican principle. Other authors have since examined 
this issue in the framework of the resolution of the dark energy problem with zero-A LT 
models (see references in, e.g., the review by Marra and Notari (2011)). 
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Quartin and Amendola (2010) performed a detailed comparison of three special 
LT void models with the redshift drift that could be measured in the future by the 
CODEX experiment, which is a high-precision spectrograph proposed for the European 
Extremely Large Telescope. They found that, for a mission duration of 15 years and 
for sources with redshifts 3 < z < 4, the drift could be around —2 x 10~^, which they 
claim might be measurable by this kind of very high-precision device. 

All these authors thus propose that the LT models might be put to the test by such 
an experiment. This is correct in principle. However, we recall here that in Sec. [T]we 
stressed that these models are not exact representations of our Universe. Even in the 
radial direction, the inhomogeneities are somehow smoothed on observable scales. This 
means that these models cannot he considered as being capable of exactly reproducing 
the features of the observed Universe on all scales, and, in particular, at the scales of 
the local Universe, since the observations show that this local Universe is not spherically 
symmetric around us as would be true if the models were exact. The very small values 
calculated for the redshift drift over a timescale of a decade or so might be used to 
put the different models to the test if these models were exact representations of the 
Universe in the region of the light cone located between us and the sources. Since this 
is not the case, either for LT models or FLRW ones, we suspect that the thickness of 
the lightcone, i.e., the time elapsed between two measurements of the redshifts, needed 
to properly test these models will be actually longer. 

Moreover, some large hump LT models have been shown to be able to exhibit a 
positive redshift drift as ACDM universes do (Yoo et al. 2011). This result sets a new 
challenge for the forthcoming redshift drift experiments, which will have to be more 
precise than expected to distinguish between inhomogeneous and homogeneous models. 
The mere sign of the effect will indeed not be enough to draw any strong conclusions, 
but its amplitude will have to be measured with sufficient accuracy. 

5. Conclusion 

Among the proposals put forward to explain the observed SNIa luminosity without 
resorting to dark energy, a number of exact inhomogeneous models have been proposed. 
The simplest and mostly studied ones belong to the Lemaitre-Tolman class with a 
central observer. We have once more stressed in this paper that this class must not 
be considered at face value but as a first step towards developing more sophisticated 
models. Their use must indeed be viewed as a mere mathematical simplification. The 
energy density is smoothed out over angles around us (Celerier et al. 2010, Bolejko and 
Sussman 2011, Bolejko et al. 2011), i.e. only the radial inhomogeneities are taken into 
account. This is, of course, a very rough approximation which has and will be improved 
in the future using more sophisticated models to deal with the issue. 

Nature does not create objects that fulfill mathematical assumptions with perfect 
precision. Objects in mechanics or electrodynamics that are described as spherically 
symmetric have this symmetry only up to some degree of approximation. An "ideal gas" 
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in thermodynamics is nearly ideal only at sufficiently low pressure. An "incompressible 
fluid" and so on. However, physicists have always employed simplifying ansatz 
to be able to complete their calculations, as in the case of exact models that are 
spherically symmetric on all scales. It is therefore misleading to test these models 
using methods that are mainly designed to test their simplifying assumptions such as 
spherical symmetry or the central spatial location of the observer (a violation of the 
so-called Copernican principle, which, as we have stressed in Sec. [H is only Copernican 
in space, not in the full spacetime). 

The same remark can be made about FLRW models. They are valid and their 
robustness has been verified on very large scales where both the curvature and the density 
contrasts remain small. However, they no longer provide a valid way to represent our 
Universe on the scales where these contrasts become too large to allow a perturbative 
method to converge. We claim that this is the case in the spacetime regions where the 
supernovae are observed. 

We hope that we have succeeded in clarifying these issues and explaining why some 
statements that can be found in the literature are wrong or vague. 

In Sec. [31 we have shown analytically that the density profile of a LT model 
reproducing the luminosity distance-redshift relation of the ACDM model is a local hill 
in the flat case E = Q and a local void in the simultaneous bang case, where = const. 
This strengthens the results obtained by Celerier et al. (2010), where the model is a 
mixture of growing and decaying modes reproducing a few features of the ACDM model 
and where the density profile is a hump near the observer. Moreover, when compared to 
the results of Iguchi et al. (2002), this shows that the energy density functions exhibit 
the same kinds of profiles for these pure cases in redshift space and on the spacelike 
hypersurface t = t^. 

A few tests applied in the literature to the LT models, and chosen among the most 
popular ones aiming at ruling out these models, have been discussed in Sec. HI 

We have shown that the kSZ effect is irrelevant to put the LT approximation to the 
test. Since, when using these models, we do not claim we are physically at the centre of 
any spherically symmetric universe but that we have merely smoothed out the angular 
inhomogeneities around us, the observer in the distant source is also not physically off- 
centre. Therefore, if she uses the same LT approximation as ourselves, she should not 
observe a large dipole in the CMB and the lack of kSZ effect seen in the telescopes 
is explained. The same remark applies to tests based on the cosmic parallax and the 
spectral distortion of the CMB. If the measurement of the cosmic parallax can provide 
an interesting evaluation of the general anisotropy of the Universe, it can say nothing 
about LT approximations. What these three tests can tell us is merely what we already 
know: that the Universe is not spherically symmetric around us. Since LT models are 
models of the Universe that are isotropic by construction, hence not exact, it is normal 
that they should be ruled out by such tests merely designed to test spherical symmetry 
or anisotropy. In contrast, it would indeed be unusual if they were to pass these tests. 

Conversely, the measured value of Hq provides a good test of the models. However, 
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in contrast to what has been advocated by Riess et al. (2011), and widely behaved 
throughout the community of cosmologists who are not experts in inhomogeneous 
models, there are LT models, even void ones, that are compatible with these values. 
Some have been ruled out, but not all. 

The most robust way to test cosmological models would be the measurement of the 
redshift drift. This could, in principle, help us to distinguish between the models since 
both the sign and amplitude of the drift depend on the model considered. However, we 
suspect that this test will need many more years of experiment before it can be said to 
be conclusive than currently advocated in the literature. 

Another criticism of the LT models, which is philosophical and not scientific, 
cannot be applied to their followers, the Swiss-cheese or meatball models, since they 
actually comply with the so-called Copernican principle. Moreover, these models tend 
to represent more physically the observed Universe with its voids and structures. The 
first Swiss-cheese models made available in the literature use the LT solution to represent 
the vacuoles spread in the background. They succeed in reproducing the luminosity- 
distance-redshift relation of the ACDM model, provided that the size of the holes is some 
hundred Mpc large. Another attempt was made by Bolejko and Celerier (2010) with 
Szekeres vacuoles where axial symmetry was used as a first simplifying assumption. In 
the future, more general Szekeres Swiss-cheese or meatball models should be constructed 
to deal with the issue. 

However, all known exact solutions of the Einstein equations that can be of 
cosmological use possess some symmetries or quasi-symmetries. The only way to 
overcome these shortcomings is to obtain a fully operational, exact, and inhomogeneous 
solution of these equations. This can only be achieved using numerical relativity and we 
suspect that this will be the new way of dealing with cosmology in the years to come. 
It will be therefore relevant to confront the models provided by these methods with 
well-designed cosmological tests. 

However, if, at the end of the day, A appeared to be nonzero or a dark energy 
component should come into play, such a result, to be robust, should be obtained with 
methods well-adapted to take the inhomogeneities into account. These could be, either 
a well-designed averaging scheme, or the use of adapted exact inhomogeneous solutions 
of General Relativity, probably numerically computed. 
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